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This paper describes numerical implementation of a newly developed hybrid model, T-blob/T-TAB , into an 
existing computational fluid dynamics (CFD) program for primary and secondary breakup simulation of 
liquid jet atomization. This model extend two widely used models, the Kelvin-Helmholtz (KH) instability of 
Reitz (blob model) and the Taylor-Analogy-Breakup (TAB) secondary droplet breakup by O’Rourke and 
Amsden to include turbulence effects. In the primary breakup model, the level of the turbulence effect on the 
liquid breakup depends on the characteristic scales and the initial flow conditions. For the secondary 
breakup, an additional turbulence force acted on parent drops is modeled and integrated into the TAB 
governing equation. Several assessment studies are presented and the results indicate that the existing KH 
and TAB models tend to under-predict the product drop size and spray angle, while the current model 
provides superior results when compared with the measured data. 


Nomenclature 

a Radius of blob or parent drop 

Bo Constant (0.61) 

B| Constant (10.0) 

C, Constant ( B 0 /3.726B, ) 

C b Constant (1/2) 

Co Discharge coefficient of injection nozzle; constant used in equation (9) 

Cp Constant (1/3) 

Cjet Constant (10) 

C k Constant (8) 

C, Empirical constant involving turbulence force 

c t Weighting parameter associated with turbulence motion 

c w Weighting parameter associated with surface motion 

D Diameter of injection nozzle 

E Term associated with energy 

k Turbulence kinetic energy per unit mass in the k — £ model 

k, Turbulent kinetic energy 
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IQ Loss coefficient due to nozzle inlet geometry 

K £ Turbulence constant (0.27) 

L Length of injection nozzle 

Ly, Characteristic length scale associated with surface motion 

Lt Characteristic length scale associated with turbulence 

r p Radius of a parent drop used in the secondary breakup formulations 

Radius of a product drop used in the primary breakup formulations 
r t Radial length scale associated with turbulence motion 

r w Radial length scale associated with surface wave motion 

s Contraction area ratio of injection nozzle 

SMR Sauter-mean radius of product drop 

t Time 

U Liquid jet injection velocity 

W Relative drop velocity with respect to local gas velocity 

y Non-dimensional parameter of distortion displacement (£/C b r p ) 

Greek Letters 


A Wave length of fastest growing wave 

p Density 

G Surface tension coefficient 

£ Dissipation rate of turbulence kinetic energy per unit mass 

T Time associated with surface wave motion 

T w Characteristic time scale associated with surface wave motion ( a/ ) 


T t Characteristic time scale associated with turbulence 

\X Molecular viscosity 

CO Angular frequency 

£ Radial cross-section change of drop from its equilibrium position 

Q, Maximum wave growth rate 

Subscripts 

n Initial value 


g Value of ambient gas 

i Value of liquid phase 

t Parameter associated with turbulence 

vv Parameter associated with surface wave motion 

Superscripts 
o Initial value 

. First derivative with respect to time 

.. Secondary derivative with respect to time 


I. INTRODUCTION 


A TOMIZATION of liquid jet into a gaseous surrounding plays an important role in combustion dynamics of 
many industrial and propulsion systems. Jet inertia, aerodynamic forces, as well as the surface tension all 
contribute to the liquid stream breakup process. In addition, the turbulence behavior inside the liquid jet also 
contributes to the jet disintegration. Often the geometrical sharpness of the injection nozzle inlet, along with 
appropriate flow conditions, can create cavitations inside the nozzle. The collapse of this cavitation can generate a 
flow fluctuation, leading to a more aggressive disintegration of the liquid jet. Traditional spray atomization models 
such as the KH instability model of Reitz [1], the TAB of O’Rourke and Amsden [2], and some enhanced models 
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[3-4], did not account for the liquid turbulence effects. Recently, Nishimura et al [5] developed a 
phenomenological cavitation model in which the liquid turbulence motion within the injector was formed by the 
bubble collapse and the fluid turbulence motion. Their model was applied to the primary breakup regime. 
Considering the nozzle exit turbulence conditions of diesel sprays, Huh et al [6] proposed a scheme taking into 
account two independent mechanisms, wave growth and turbulence in the atomization process. The turbulence is 
characterized partially by the injection nozzle geometry while the wave growth is derived from the KH instability 
theory. The aforementioned models were formulated with semi-empiricism. Several efforts of using sophisticated 
numerical approaches have been performed in modeling the detailed turbulent flow fields in the liquid and gas 
during the atomization process. In the Large Eddy Simulation (LES) spray modeling of Apte et al. [7], a stochastic 
approach for droplet breakup taking into a range of product-droplet sizes was developed for modeling of secondary 
breakup regime. De Villiers et al [8] applied LES to resolve the primary atomization under more realistic operating 
conditions found in Diesel engines. These methods have a real potential of providing a complete physical 
description of the liquid jet breakup with minimum assumptions; however, they require submicron spatial elements 
in size and pico-second in time steps to properly predict the atomizing sprays at the high-velocity injection 
conditions. Consequently, grid mesh size and the considered physical domain must be taken care of so that the 
computational time and memory storage requirements can be manageable. At the present time they are still too 
expensive and generally impractical in terms of computational time and power requirements for engineering 
calculation applications. Hence, the engineering analysis and design of the liquid spray devices still must rely on 
phenomenological engineering models. 

Based on two widely used models, Reitz’s primary atomization (blob) [1] and the Taylor- Analogy-Break (TAB) 
[2] secondary droplet breakup, a hybrid T-blob/T-TAB model was recently developed by Trinh and Chen [9] to 
include turbulence effects in both primary and secondary breakup regimes. The “hybrid” model is used to explicitly 
differentiate between the primary and secondary breakup mechanisms that occur in different regimes. The 
approaches described in [5-6], and others [10-11] are hybrid models. Here, in the primary breakup model, T-blob, 
the level of the turbulence effect on the liquid breakup depends on the characteristic time scales and the initial flow 
conditions. This treatment offers a balance of contributions of liquid turbulence and perturbed surface wave physical 
phenomena on the liquid breakup process. For the secondary breakup, T-TAB, an additional turbulence force acted 
on parent drops is modeled and integrated into the TAB governing equation. The drop size formed from this breakup 
regime is estimated based on the energy balance before and after the breakup occurrence. The turbulence energy is 
also considered in this process. Detailed model development and theoretical formulations, as well as several 
fundamental test cases and energy budget calculations were described in [9] for simplified flow fields to validate 
several model coefficients. The purpose of this paper is to incorporate this newly developed hybrid 
phenomenological T-blob/T-TAB model into a Computational Fluid Dynamics (CFD) methodology based on 
Eulerian (gas)-Lagrangian (liquid) formulation for complicated injector liquid-atomization and spray simulation. 


II. MODEL DESCRITION AND IMPLEMENTATION 

Detailed development and theoretical foundation of T-blob/T-TAB model was described in [9]. Only the final 
forms of the model are summarized here. 

T-blob Model 

In this model, the drop breakup rate and the mean drop size are given by: 


da 

~ck 
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r p <a. 


( 1 ) 


r P = 


(l-c t )r t +c t r w 


( 2 ) 


3 

American Institute of Aeronautics and Astronautics 



The characteristic length scale L* and time scale T w are associated with the droplet surface wave instability. These 
two scales along T and C a are formulated from the “blob” model [1] as: L w = A , x w = a/AQ , x = 3.726B,a/A£2 

The turbulence characteristic length scale Lt and time scales T t were derived using the 


and C = — — — 

1 3.726B, 

analytical solution of the k - £ turbulence model: 
X, =<+0.0828t 


L. =L® 


1 + 


0.0828t 
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\ 0.457 


(3) 

(4) 


The time t is counted from the time at which the parent drop leaves the injection nozzle exit. The initial turbulence 
length scale L° t and time scale X° t are evaluated from the initial turbulent kinetic energy k° and its corresponding 

dissipation rate E° t at the injector exit: 


S n 


T ? = c ti “o' ’ where c „ = 0.09 • 


The initial turbulent kinetic energy and its corresponding dissipation rate 8° can be approximated as 


‘ 8L/D 

= K — 
1 2L 
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LW 




(5) 


(6) 


with ^=0.27- The velocity U is the liquid velocity at the injection nozzle, which has the length, L, and the 

diameter, D. The discharge coefficient, the loss coefficient due to the nozzle entrance sharpness, and the 
downstream-to-upstream contraction area ratio of the injection nozzle are represented by C d , Kc, and s, respectively. 
The parent drop would no longer strip its mass to create the product drop when the parent drop radius a is less than 

the radius r p of the product droplet. In Equation (2), determination of the radius r w associated with the wave motion. 


follows the original formula of [1]; and the value of r t , as well as weighting coefficients can be found in [9]. Along 

with the inclusion of the turbulence effect on the primary atomization process, the subject phenomenon also is 
transitioned to the secondary droplet breakup model, the T-TAB model. The required initial velocity fluctuation 
quantity of the product drops right after their formation was obtained by examining the energy conservation during 
the primary breakup process. By equating the change in the energies of the parent drop with the energy of the 
product drops, the turbulent kinetic energy for the product drop was derived [9] and was used to estimate the initial 
velocity fluctuation tor the secondary breakup process, which is discussed next. 


T-TAB Model 

To account for turbulence effect, the original TAB model was modified to: 


C b p, r p 2 ‘CftLk, 


-|C. 


y'-C k 


P. r P 


ycAi 
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(7) 


Except for the second term on the right hand side, equation (7) is the same as the governing equation of the TAB 
model. Again, reference [9] should be consulted for detailed derivation of this model. To estimate the post-breakup 
drop size, a similar method as used in the TAB model is employed for the calculation of Sauter mean radius (SMR): 

4jtr„ 3 (T ra\ 

SMR • w 
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Based on the conservation principle, the initial turbulent kinetic energy ( k 0 ) and dissipation rate ( £ 0 ), in Equation 

(7), of the new drop were formulated [9]. in Equation (8) is the total energy of the original drop prior to the 
breakup and is determined from the original TAB model. 

T-blob Model Implementation 


The liquid is assumedly injected in a form of “blob” parcels containing spherical drops with their initial size 
equal to the nozzle exit dimension. The computational approach for the primary atomization proposed by Reitz [1] 
is adopted in the present study. Extra relations representing the turbulence are described here. The initial turbulence 
quantities are estimated based on Equations (5) and (6). Values of the turbulence characteristic scales are then 
calculated from these turbulence values. Individual parcel is tracked at each time step, and change in “blob” drop 
size is also determined using Equation (1). Next, the mass stripped from the parent drops is computed and is formed 
into a new product drop parcel. The primary breakup of the blob drops continues until their size reduces to the 
dimension of the product drops. 

It is not necessary to create a new parcel in each computational time step. For the purpose of saving computer 
memory and computational time, the stripped mass would be accumulated until certain criteria are met. Then, a new 
parcel is added to the computations and its drop size is calculated from Equation (2). The criterion of forming a new 
parcel varies among authors. For instance, Reitz [i] defined the creation of a new parcel when the accumulated 
mass had exceeded 3% of the original parcel mass; while Huh et al. [6] used the value of 10%. It was found that at 
least for the test cases of interest, the computational results from theses two criteria are not much different, although 
the number of parcels in the computational domain varies considerably. In CFD-ACE+ [13], a new parcel is created 
when the number of product drops is greater than 20% of the ones in the parent. In the present study, the initial 
turbulent kinetic energy of the product drops is estimated from the energy conservation formulation described in [9]. 
When this estimated energy is less than 0.1% of the kinetic energy of the droplet motion, it is assumed that the 
turbulence motion within the droplet becomes insignificant. Then, no turbulence is considered in the secondary 
breakup for subject droplet, and the classical TAB model is employed. Otherwise, the drops would be considered 
for the secondary breakup with the turbulence effects that are described in the next subsection. 

T-TAB Model Implementation 

While the parcels produced from the primary atomization are continuously tracked in the Lagrangian coordinate 
system, they are also considered for the secondary breakup, which is governed by Equation (7). This equation is 
solved in the same way as the classical TAB model [2] but using the finite difference method since its exact 
analytical solution does not exist. 


To develop the finite difference formulation for Equation (7), it is convenient to define a new set of coefficients 
for this equation as follows: 


C,= 


C-b Pi r p 


c, =c. 


' 2 cv 
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By substituting the above coefficients, Equation (8) becomes 

y = C,+C 2 y ‘-C 3 y-C 4 y 


In contrast to the TAB model, an exact solution can not be obtained for since it contains an extra non-linear term 
(second term on the RHS). Instead, a second order finite central difference scheme in a time step was used [12]. The 
deformation displacement y is computed at every time step. The secondary breakup assumedly occurs when y 
approaches 1. Then the product drop size is calculated from Equation (8). At the same time, its initial turbulence 
quantities are estimated [9]. Again, the new product drops are continuously tracked for possible future secondary 
atomization until they leave the computational domain. When the product drop size is smaller than the Kolmogorov 
length scale, the turbulence activity of the new drop is no longer considered. In this case, the present model would 
recover the classical TAB model. 


The T-blob/T-TAB model was incorporated into the CFD-ACE+ [13] program for two-way gas-liquid 
atomization calculations. This code, together with its liquid spray module, solves the equations of steady-state as 
well as transient flow with chemically reactive fluid dynamics and droplet spray. The governing equations and the 
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numerical solution method are discussed in [13] and are only described briefly here. The Navier-Stokes equations 
and RNG (Renormalized Group) k — € turbulence model for a gas phase are solved with a finite volume method. 
A time-marching with a high-order accurate explicit scheme is employed in the transient calculation. With an 
assumption of dilute spray, parcels of drops are tracked in a Lagrangian framework. The parcel represents an 
ensemble of identical drops occupied in the same space at a given time. Subsequently, the spray properties at each 
point are described by statistically sampling the spray parcels. The drop parcels exchange mass, momentum, and 
energy with the gas through source terms in the Navier-Stokes equations. 


III. RESULTS and DISCUSSION 

Numerical simulations of several test cases as listed in Table 1 are performed and compared with available 
measured data. Variations in the predictions from the two classical KH and TAB models and the present models are 
also examined. The two-dimensional axisymmetric flow of a cylindrical liquid jet located at the center line of the 
computational domain is simulated. It should be noted that the current version of CFD-ACE+ does not have 
collision/coalescence models. Due to the observation of these phenomena [15], a module was written [12] to 
compute the droplet collision and coalescence based on model of [14]. The time-transient liquid jet breakup and the 
ambient gas behavior are computed for an incremental time step, At, of 2.5xl0' 6 sec. The spray trajectory is tracked 
in the Lagrangian coordinate system. A different time step, which is normally smaller, is applied for this coordinate 
system. CFD-ACE+ recalculates this time step continuously based on criteria of pre-assigned maximum allowable 
changes in drop velocity and diameter. The Lagrangian tracking time step used in the simulations for this study 
ranges from 1/3 to 1/10 of At. It should be noted that a solution sensitivity study for various numerical parameters, 
such as grid mesh size, time step values, etc. has been performed [12] prior to obtaining the results of this 
investigation. 


Table 1. Test cases and measured data used in the computations 


Case 

H-l 

H-2 

H-3 

K 

Nozzle Diameter (mm) 


0.3 


0.24 

Ambient Gas 


Nitrogen 


Nitrogen 

Ambient Pressure (MPa) 

1.1 

3.0 

5.0 

2.17 

Ambient Temperature (K) 


298 


298 

Density (kg/m 3 ) 

12.36 

33.70 

56.17 

24.51 

Liquid Fuel 


Diesel Fuel 


Diesel Fuel 

Density (kg/m 3 ) 


840 


840 

Viscosity (Kg/m’ s) 


2.9xl0' 3 


5.0xl0 3 

Surface Tension (N/m) 


2.05x1 O' 2 


2.06x1 0‘ 2 

Injection Velocity (m/s) 

102.0 

90.3 

86.41 

133.81 

Initial Turbulence Quantity 





Kinetic Energy (m 2 /s 2 ) 

2.88x1 0 2 

2.26xl0 2 

3.94x10* 

5.22X10 2 

Dissipation Rate of Kinetic Energy 
(m 2 /s 3 ) 

1.06x10® 

7.34xl0 7 

1.22xl0 7 

3.14x10® 

Reference 

Hiroyasu et al. [15] 

Koo [16] 


The spray tip penetrations from the exit port into the ambient gas for test cases H-l and H-2 are plotted in Figure 
1 as a function of time. Both predictions from the T-blob/T-TAB and KH/TAB are similar. Nearly straight lines of 
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the tip penetration curves at the beginning indicate that the tip velocity of the liquid jets remains almost constant. 
Once the jet is fully converted into droplets, the surface area of the spray nose is enlarged, and its speed is 
subsequently reduced due to greater aerodynamic force and drag. Hence, the penetration curve has a smaller slope 
at the later time. Furthermore, the results also reveal that the spray penetrates at faster rate into a lower pressure gas. 
In comparison with measurements, the tip penetrations in both cases are slightly under-predicted with analytical 
models at the beginning of the spray; however, they are over-predicted somewhat at the later time. Overall, the 
computational results show a close agreement with the experimental data. Variations in the spray angle due to 
different ambient gas pressures for test cases H-l, H-2, and H-3 are depicted in Figure 2. This angle is a global 
parameter used to characterize the radial expansion of the spray. For an axisymmetric case this angle, as shown in 
Figure 2, is formed by the two opposite tangent lines drawn along the spray envelop. The results portray higher 




Figure 1. Comparison of tip penetration. Figure 2. Variations in spray angles vs. back P 


back pressure with a larger radial expansion of the spray. A further examination of the atomization modeling 
suggests that the radial product drop velocity is related to the wave length and its maximum growth rate. The values 
of these parameters, however, are dependent on the liquid-to-gas density ratio. Consequently, higher gas pressure 
would result in a larger radial-to-axial velocity ratio and a wider spray angle. In comparison with the results of the 
KH/TAB models, a favorable larger spray angle is predicted by the present models. As stated in [9], the T-blob 
model estimates considerably larger product drop sizes than the KH model, which leads to the prediction of a 
stronger radial momentum. Consequently, liquid drops are spread outward more in the T-blob model when 
compared to the predictions provided by the KH models. As shown in Figure 2, a better agreement with the 
measured data is obtained with the present models. 

The computed spray shapes and tip penetrations from the present models are displayed in Figure 3 for three 
different back pressure conditions after 2.5ms of liquid injection. It should be noted that the injection velocity is 
higher for a lower back pressure condition. Hence, the spray tip penetrates faster into a lower pressure environment, 
although the drag increases on the periphery of the spray contour. This phenomenon leads to parcels containing 
small drops being dragged into the wake behind the tip region, which clearly can be observed in Figure 3 (a) and (b). 
It is also interesting to note that larger drops appear on the spray tip as well as on the peripheral regions near the tip. 
Generally, large drops carry a stronger momentum, which decays less than their smaller drop counterparts. Hence, 
the large drops tend to travel faster downstream than the smaller drops do. In addition, large drops found in these 
regions may reflect the effect of drop coalescences. The coalescence is more pronounced in an elevated gas 
pressure. As illustrated in Figure 3(c), a larger drop size can be observed in the tip region at higher ambient 
pressures. 

The following results are collected from numerical simulations of test case K (see Table 1). This test case is 
selected to evaluate the atomization models because of its available measured drop size data. In the experiment, 
Koo [16] measured diesel sprays with a Phase/Doppler Particle Analyzer (PDPA) system. In general, a PDPA can 
recognize only a defined range of droplet diameters due to digital processing hardware limitations. Normally, the 
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ratio of the largest observable droplet to the smallest observable droplet cannot exceed 35. Regarding to this 
experiment, the instrument was set to measure drop sizes from 1.8 to 246 pm. It should be noted that the liquid fuel 
injection was driven by a cyclic pumping system for which the mass injection rate is no longer constant. To 
simulate more realistically the inlet flow conditions, the transient injection velocity profile reported by Koo [16], as 
shown in Figure 4, is utilized in the computations. 
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Figure 3. Variations in spray shape, tip penetration, and drop size at t=2.5ms due to different back pressures, 

(a) 1.1 MPa; (b) 3 MPa; (c) 5 MPa 



Similar to the previous computational cases, this test case is simulated also for the transient flow field conditions 
with At=2.5xl0' 6 s. It should be noted that the present simulation utilizes the RANS approach to obtain the mean 
velocity field of the gas phase. For statistically non-stationary flows, which were investigated experimentally by 
Koo [16] in this case, the interpretation of the averaging process is the ensemble averaging; thus the RANS 
formulation is justified for the statistically transient case as indicated by the injector inlet condition of Figure 4. The 
predicted spray tip penetration and the measured data are presented in Figure 5. Also shown in this figure are the 
grid-refinement comparative results. With finer grids packed in the near nozzle region, the spray seemed to be 
“denser” and penetrate further into the gaseous environment. This is the general shortcoming of Eulerian-Lagrangian 
approach utilizing particulate “point source” assumption within the Eulerian grids [17]. Extremely fine grids can not 
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be used. Within certain grid density, the spray exhibits similar behavior. The predictions of this test case, made with 
both the KH/TAB and present models, are slightly lower than the measurement at the beginning of the injection; 
however, the predicted curves approach the experimental data at a later time. It is noted that Koo reported the spray 
penetration length from the beginning of the injection to the time of 1.5ms; while the computation continues to 
3.0ms. In general, both the predicted tip penetrations are similar, and their values are in reasonable agreement with 
the experimental data. Figure 6 displays the simulated spray shapes along with the photographs recorded at several 
time intervals in the experiment. In comparison with the photographs, the simulated spray penetrates with a slightly 
shorter distance into the ambient gas than the actual distance depicted from the photographs initially; however, the 



Figure 5. Predicted and measured spray tip penetration for test case K 



Figure 6. Simulations and photographs of case K; (a) t=0.203 ms; (b) 0.601 ms;(c) 1.205 ms;(d) 1.8 ms 
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Figure 7. Drop size distributions from the predictions and measurements for test case K 
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Figure 8. Drop size distributions from the predictions and measurements for test case K 
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predictions improve significantly afterwards. These results also reflect the same assessment of the tip penetration 
comparison shown in the previous Figure 5. It is also interesting to observe the degree of the spray being spread 
laterally and its tip shape. All the models predict a similar radial expansion of the spray for all considered time 
intervals, and they are compared reasonably well with their photographic counterparts. However, the numerical 
simulations display the tip profile more like a conical shape with parcels containing large drops rather than a round 
tip shown in the photographs. Since the momentum of these large drops is degrading much less than the one of the 
small drops, the large drops move faster downstream than the smaller drops. Therefore, the moving motion of the 
large drops creates the appearance of the conical shape at the tip. A series of the computational cases has been 
conducted to examine the tip shape variation for different initial spray angles. The analysis suggests that the spray 
tip head-end changes to a round-shape and does not penetrate into the ambient gas as deep when increasing the 
initial spray angle. The tip variation analysis is not addressed in detail here, since it is not a main focus in this study. 

Further detailed drop size comparisons between the predictions and the measurements are presented in Figures 7 
and 8. Figure 7 reveals the drop sizes recorded at three axial locations along the center line of the spray, which are 
identified in Figure 7(a). The results indicate that the jet, with an initial diameter of 240 pm, breaks up into 
relatively large drops along the center line. Their product drop diameters range from 100 to 200 pm at the stations 
of x=10 and 20mm (Figures 7(b) and (c)). These results suggest that the primary breakup process is dominant in 
these regions. Further downstream, the drops undergo the secondary atomization which generates smaller drops. 
Figure 7(d) exhibits drops of 20 to 50 pm in diameter at the x=60mm station. It should also be noted that the fuel 
injection cycle for this test case is complete at 1.77ms; therefore, the drop sizes shown in Figures 7(b) and (c) also 
reflect this end period at which the small drops are registered. In Figure 7(b), the effect of finer grid within this 
region is also observed. As discussed earlier, the finer grid system seems to see “more” droplets even the spray has 
been terminated (after 2 ms). The small drops recorded after 1.77ms are formed because of the small injection rate at 
the tail end of the injection cycle, as well as the small drops lagging behind in the spray motion. In general, the 
predictions from the present models are in good agreement with the measured data. In contrast, the HK and TAB 
models under-predict the drop size for stations of x=10 and 20mm; however, their predictions are improved 
significantly at the x=60mm station. This observation is consistent with earlier results in [9]. The results also show 
that the KH model predicts generally smaller drop sizes than the T-blob model. Figure 8 displays the drop sizes at 
three off-centerline locations [x, y] shown in Figure 8(a). For the location at [10, 0.5], the measurement registered 
relatively large drops, the diameters of which are ranging from 30 to 180 pm with a random order in time. Similar 
to the results shown in the previous figure, the smaller drops are also observed after the end of the injection cycle at 
t= 1.77ms. The present model predictions capture large drops intermittently at certain short time intervals; however, 
both the sets of the KH/TAB and T-blob/T-TAB atomization models project typically a much smaller drop size than 
the measurement portrayed. In other words, the CFD simulation indicates that the secondary droplet breakup 
process takes place at this location but not the primary breakup. The drop diameter at [10, 1.5] and [30, 2] in Figures 
8(c) and 8(d), respectively, varies from 20 to 50 pm. This suggests that considerably smaller drops appear on the 
spray periphery and their presence in this outermost region forms the shape of the spray. Furthermore, since the 
drop sizes at these locations are relatively small, the secondary droplet breakup occurs mostly in these regions. This 
assessment is consistent with the observation of Koo in his measured data analysis. He concluded that the small 
droplets in the periphery of the spray are not those directly injected from the nozzle. The small droplets from the 
nozzle travel much slower than the large droplets, which are produced in the primary breakup; hence, they can not 
maintain the spray shape observed in the experiments without the appearance of drops formed by the secondary 
breakup mechanism. It is interesting to note that the simulation with the KH and TAB models has a limited number 
of parcels and no parcel tracked at [30, 2] and [10, 1.5], as shown in Figures 8(c) and 8(d), respectively. Again, this 
examination reconfirms the display in Figure 2 that the spray angle predicted by the classical KH and TAB models 
is a little smaller than the actual angle. In the present model, the numerical simulation captures droplets at [10, 1.5] 
and [30, 2] only in a short time span, and the predicted drop sizes at this interval are close to the measured data. 

IV- CONCLUSIONS 

We have successfully implemented a newly developed hybrid T-blob/T-TAB atomization model into an existing 
CDA-ACE+ code for simulation of liquid jet atomization using two-way Eulerian-Lagrangian method. A drop 
collision and coalescence model [14] has also been incorporated into the ACE+ code to provide more realistic 
simulations. Assessments of the new model have been performed against two experimental studies and several 
observations and conclusions can be drawn as follows: 
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1. In general, the CFD simulations are able to capture similar spray structures observed in the experiments. 
The spray tip penetration lengths as well as the appearances of large drops at the spray tips, which are 
predicted from the two sets of the KH/TAB and T-blob/T-TAB atomization models, are comparable to the 
measured data. However, the KH and TAB models estimate smaller spray angles for the test cases H-l, H- 
2, and H-3 than the measurements indicate. 

2. For the test case K, the drop measurements and the predictions indicate that small drops created from the 
secondary breakup appear in the outermost regions of the spray. Hence, this breakup mechanism is 
responsible for the radial expansion of the spray. These drops then form the profile of the spray. On the 
other hand, the drops on center line of the spray are relatively large. Furthermore, the appearance of large 
drops at the spray tips may result from drop coalescence. In addition, the strong momentum of the large 
drops allows them to travel ahead downstream. 

3. In summary, for test case K, the CFD simulation using the HK and TAB models tends to under-predict the 
drop size of the spray, while it is able to capture reasonably the liquid spray structure. On the other hand, 
the present models predict comparable drop sizes, as seen in the measured data. Overall, the predicted 
results from the present models do agree reasonably with the experimental data. 
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